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Abstract 

In this paper we investigate the subleading twist corrections to the photon energy 
spectrum in the decay Aj, — » X s ~/ and to the lepton energy spectrum in A& — » X u £i>£. 
As a first step we rederive the matching coefficients of the subleading twist us- 
ing a much simpler method. Parametrizing the matrix elements for A& and using 
reparametrization invariance, we show that the energy spectra close to the endpoints 
are given in terms of only a single universal function. 



1 Introduction 



The existence of baryons with bottom quarks is well established by now, and at future 
facilities a detailed study of the decay modes of such states will become feasible. However, 
the B factories are running below the threshold for the production of these modes, which 
means that the measurements we are going to discuss will be - at least in the foreseeable 
future - the domain of hadron colliders. 

From the theoretical point of view, the ground state bottom baryon A& is the simplest 
heavy state, although it contains three quarks. The reason is that the spin of the heavy b 
quark decouples in the limit m;, — ► oo and thus the light degrees of freedom have to be in 
a spin- zero state pQ . In other words, the spin of a A& is the spin of the heavy quark, which 
in principle allows interesting polarization studies. Some of these studies require polarized 
Ab's which may become available at GIGA-Z j2], a high-luminosity linear collider running 
at the Zq resonance. 

In the present paper we investigate energy spectra of inclusive heavy-to-light decays 
of the A ft , in particular the photon spectrum of the rare decay A b — > X s ^ and the lepton 
energy spectrum of the semileptonic decay A5 — > X,J,v^. The situation here is very similar 
to the one in the decays B — > X s ^ and B — > X u £uf. The photon spectrum in B — > X s 7 
is mainly concentrated in the endpoint of small hadronic mass of X s , where the strict 
1/m expansion breaks down. This forces us to switch to an expansion in twists similar 
to deep inelastic scattering. The leading twist contribution is given in terms of a light- 
cone-distribution function which can be formally written as a forward matrix element of a 
non-local operator This function describes all inclusive heavy to light transitions and 
is universal. Applying this to b — > utvi and b — > 57 one may obtain a model independent 
determination of | / Vt s | 

Recently subleading terms have been included into the analysis of B meson decays [5 , 
and the purpose of the present paper is to include subleading terms for the inclusive decays 
of heavy baryons. The main result of the paper is that only a single distribution function 
is needed to parametrize the differential decay rates to subleading order; in other words, 
the endpoints of the photon energy spectrum and the charged lepton energy spectrum in 
At — ► Xsj and Ab — > X u ivt are given in terms of one universal function up to terms of 
order l/m^. This has to be compared to the case of the corresponding B decays, where 
three universal functions are needed [5]. 

In the next section we consider the matching including the terms of subleading twist. 
We reproduce the results of [5] using a much simpler method based on a background field, 
in which the light quark propagates. In section 3 we parametrize the forward matrix 
elements for the non-local operators between Ab states. In section 4 we discuss our result 
and conclude. 
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2 Matching on non-local operators to order 1/m 



The tree-level matching of full QCD to the non-local operators to order l/m b has been 
performed for both b — > s r y jH] and b — > utvi [HI E] and the results of these papers can be 
used also for the case of the Ab. However, we find it useful to rederive the results of [5] in 
a different way, which is very transparent and is manifestly gauge invariant. The method 
resembles the one used in jH] for the shape function for mesons. 

The starting point is the usual decomposition of the heavy quark momentum according 

to 

p b = m b v + k (1) 

where the components of the residual momentum k are small compared to m b . At tree 
level, the inclusive decay b — > 57 requires to consider the propagator of the (massless) s 
quark, which is 

S = % = - (2) 

where q is the photon momentum 

In position space and in terms of the static heavy-quark field h v (x) the residual mo- 
mentum k corresponds to a covariant derivative, since the large part m b v has been removed 
by a field redefinition. This leads us to consider the following (non-local) operator 

T = h v (0)Vn -AYhJx) (3) 

where D is the (QCD) covariant derivative and T is the Dirac matrix of the corresponding 
heavy-to-light current sTb 

It has been shown in [8 that (jSJ) yields the correct tree level matching for the shape 
function and we shall use the same arguments to rederive the subleading terms. In order 
to define the power counting properly we introduce light-cone vectors as 

v = -(n + n), n 2 = = n 2 , n ■ n — 2, q — (v ■ q)n (4) 

In this way we can decompose any vector into its components according to 

1 1 ± 

SV = 2 n A^ + 2 n » n » + V ( 5 ) 

We want to expand in the kinematical region where 

(m b v -q)-n = m-n-q~ A QC d (6) 

which is the kinematical region for the shape function. In order to perform this expansion 
we write 

1 m b f + ip — ^ 



m b f) + ip — </[ {m b ^ + ii 



(7) 
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and consider the denominator, which can be written in the form 

{m b f + ip- if = m b Q + + ^{iD-, Q+} + {tD 1 ) 2 - ia^iD^iD v (8) 

where we have defined 

Q + = m b -n-q + iD + , iD + = i(n-D), iD_ = i(n-D), = g±D" (9) 

Since m b — n ■ q ~ Aqcd we have Q + ~ Aqcd and thus the term msQ+ is the leading 
term (of order m b AqcD) in the denominator. All other terms are of order Aq CD and are 
expanded as 



1 



{mj> - i + ip) 2 mQ + mQ + 2 

1 



mQ 



(10) 



+ 



mQ A 



((iD^-ia^iDJDv) 



mQ, 



+ 



Next we reinsert the numerator, which we do in a symmetrized fashion. After some simple 
algebra we get 



mi 



ip 



+ 



+ 



{mi-i + ip) 2 2Q+ 2m 2m 

1 1 w 



2m Qh 



Q+j 



(11) 



The last term can be rewritten in terms of the perpendicular components of the covariant 
derivative. We have 



1, a^ u }n^n v = 



+ - = o 



from which we derive 

{i, ia^YD,iD v = ia^n a n pl l\iD^ Q + ] + iD^iD^ 
Combining this with the other terms we get 



mf — 



71 



{mf - i + ip) 2 2Q + 2m 2m 
1 i 



4m 



W -Q- + 



2m Q A 



(12) 
(13) 

(14) 



(\g?{iDi,iDt} - \i<[iD^iD^ J- (15) 
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This result can now be used to perform the matching on the non-local operators as 
they have been defined in j^j. These operators are 



Oq{u) = h v 5(u + iD + )h v 

Of(u) = h v {iD li1 6(u + W+)}h v 

2 (u) = h v iD^, 5(u + iD+) h v 

Of (wi,w a ) = h v 6(uj 2 + iD + ){iDl, iD^diut + if> + )h v 
0^(ux, uj 2 ) = h v 5(oj 2 + iD + )[iD^_, iD±]5(ux + iD + )h v 



and 



pr 



h v 5(w + iD + )j a j 5 h v 

h v \iD^ 5(uj + iD + )}^ a 7 5 h v 

h v iDn, 5(ou + iD + ) jaTf 5 h v 

h v 5(u 2 + iD + ){iD^, iD v ± }5(u 1 + ib + )^ a ^ 5 h v 

h v 5(u 2 + iD + )[iD^,iD" ± }5(ui + iD + )~/ a "/ 5 h v 



(16) 
(17) 

(18) 

(19) 
(20) 



(21) 
(22) 

(23) 
(24) 
(25) 



The matching onto these operators is performed for the decay rate. In order to obtain 
the rate we have to take the imaginary part of a forward matrix element with the A&. 
Inserting ()15j) back into (jBJ) and taking the imaginary part we obtain for the leading term 



Im h v ——h v 



1 



7T 



/ duo -Tr [P+T^T] 5(1 - x - uo)h v 5(w + iD + )h v 

2m J \2 

~Tr [s a rfiF] 5(1 - x - uo)h v 5(w + iD+)^ b h, 

_J_ f duo[C (uo)O (uo) + C« (uj)P , a (uj)] 
2m J 



(26) 



from which we can read off the matching coefficients Cq and C§ as they have been derived 
already in [S] 1 

In the same way we obtain for the subleading contributions 



Im h, 



ME/- 1 

' 2m \ 



Wl " Q. 



(27) 



r - r7 M r - i 

Im J duoh v \^ {iD^,5(uo + iD + )}- 



-K 



x — UJ 



^^Here we made us of the decomposition formula 



P+TP+ = -P + Tr[P+T] - -s»Tr[ Slx T] = P +W P+ P+ = ~(1 + ft. 
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2m 2 



rL; ( —Tr 



P+r 7 jTj 5(1 - x - u)h v {iD^ 5(u + iD + )}h % 
s a r^T 5(1 - x - uj)h v {iD^5{uj + iD + )}^ 5 h v 
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Tr 

2 

from which we obtain the matching coefficients Cf and C^f 

A further restriction of the structure of the matching coefficients, which should hold 
even beyond tree level, is due to reparametrization invariance [9J. It turns out that 
reparametrization invariance jTTT|. ITTj relates the coefficients of Oq and 0$ v in such a way 
that only one function is needed to parametrize the effect of the two operators Oq and 0^ u . 
We shall use the relations from ^0] in the next section when we will count the number of 
of unknown functions needed to describe the spectrum of A b — » X s ^ including subleading 
contributions. 

In addition to these operators we have also to take into account the subleading terms 
of the Lagrangian. They are given as 



O t {uj) 



d x 



dt 

2tt' 
dt 



-iujt 



T[h v (0)h v (t)C 1/m {x)} 



I d " X J ^ e ~" tT M°)l^h v (t)C 1/m (x)] 



where C\/ m is the subleading contribution to the Lagrangian 



^-h v {ip) 2 h v 
2m 



(28) 
(29) 

(30) 



Once matrix elements are taken, these pieces can be interpreted as the corrections to 
the states, which means that they will appear always in a specific combination with the 
leading terms, i.e. these will not introduce any new unknown functions. 



3 Matrix Elements for Decays 

After having computed the matching we have to consider the forward matrix elements 
with a Ab. These matrix elements have to be evaluated in the static limit in which the A b 
becomes a very simple object. We obtain for the matrix elements: 

(A b (v,s)\O (uj)\A b (v,s)) = u(v,s)u(v,s)f A (u) (31) 
(A b (v, s) \P , a {u) \^b(v, s)) = u(v, s)7 a 7 5 w(f , s)f A (u) (32) 

where u(v, s) describes the static A b with velocity v and spin s. 

For the subleading operators the non-vanishing matrix elements are parametrized as 

(A b {v, s)\0?{u)\A b {v, s)) = -2u{v, s)u{v, s){v" - n»)uf K (u) (33) 
(A h (v,s)\P» a (uj)\A h (v,s)) = -2u(v,s) lal5 u(v,s)(v»-n»)ujf A (uj) (34) 
(A b (v, s)|Of (^i,^ 2 )|A 6 (t>, s)) = u(v,s)u(v,s)gfg\(u 1 ,u 2 ) (35) 
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(A b (v, s)\P^(ui,u 2 )\A b (v, s)) = u(v,s)j a j 5 u(v,s)g l ^g A (uj 1 ,uj2) (36) 
(A b (v, s)\0 T (u)\A b (v, s)) = u(v, s)u(v, s)t A (u) (37) 
(A b (v, s) |Pt,«(w) I A b (v , s)) = u{v, s)j a j 5 u(v, s)t A (u) (38) 

Consequently, only two additional functions appear at order l/m b . Furthermore, these 
subleading function are tied to the leading order function f A by reparametrization invari- 
ance, which means that it appears always the combination 

F\(u) = f A (u) + t A (u) - - 2 lMu l9A M ( (39) 

mf J \ u)\ u)% J 

with the leading order function f A . In other words, due to heavy quark spin symmetry we 
can express A b — > X s ^ and A b — > X u ti>t in terms of a single universal function, up to terms 
of order l/m^. 

Finally we discuss the moment expansion of the resulting functions. For the leading 
term we have 

/ A (u;) = %)-^>) + ... (40) 

where 

X A u(v, s)u(v, s) = (A b (v, s)\h v (iD) 2 h v \A b (v, s)) (41) 

is the kinetic energy of the b quark inside the A b baryon. For the subleading contributions 
we obtain 

2X A 

g A {uj!,uj2) = — SMSM H (42) 

t A (u) = -^V) + " " " (43) 



which reproduces the known result from [12] with the replacements Ai — > Aa and A2 — > 0: 

(44) 



|(A^.Y,0) = r„ 



„ Glc\V£V t ,\\Ctf 



with 



327T 4 

The moment expansion of the new universal function reads 



ml (45) 



F a (uj) = 5(uj) + -^5'(uj) + --- (46) 
bmi 
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4 Discussion and Conclusions 



The main result of this paper is that the semileptonic and radiative decays of a A& baryon 
into light hadrons are described by a universal function, where the corrections to this 
statement are of order \jm 2 - In this way we find the two differential decay rates: 



T Afe 
dx 



G%a\V t tV ts \\C 7 \ 2 



X u ivp 



G F | V u b 
96vr 3 



32tt 4 

2 



-777,6 



-m b 5 (2x-l)F A (l-x) 

0(1 - y - - uj)F a (u) 



(47) 
(48) 



where x = 2E 1 /n r ib is the rescaled photon energy and y = 2Ep/rrib is the rescaled lepton 
energy. Note that the first moment of F still vanishes up to terms of order l/m^ 2 , so we still 
get no contribution to the total rates. This serves a a check, since the contributions of order 
I/777 may not contribute to the rate. As for Af, — > X s j, the spectrum of the semileptonic 
decay can be expanded in terms of moments. Inserting ()46|) into the spectrum, one arrives 
at 

Aa . m , Aa 



dy 



XJPA 



0(1 -y) 



S(l-y) 



Qmh 2 



S'(l-y) + ... 



where 



Gp I V u b 1 2 



96tt 3 



-777,6 



(49) 



(50) 



With the replacements already mentioned above, this agrees with the result obtained in 



Although this may be difficult from the experimental point of view, we may investigate 
the subleading corrections to a determination of V u b from a comparison of A^ — > A s 7 
with Ab — > X u ii>i. Like in the case of B mesons, an energy cut E c on the charged lepton 
will be unavoidable to discriminate the large charm background, i.e. E c > (M 2 [A6] — 
M 2 [A C ])/ (2M[A&]). One may define observables (similar as for the decays of B mesons, see 
[T5] ) involving partially integrated rates with suitable weight functions. Since the 1/rrib 
terms appear only a kinematic factors (the pre-factors (2x — 1) in (}4Tj) and (1—u) in (|48|1) 
we have more complicated weight functions as in [13J. We define 



r u (£ c ) = 



m A /2 



dEp 



'AE l 2 -E c m h \ dT{ 



IE, 

2 

m b JE, 

from which V u b can be determined as 

Vub 



2E e 



T S (E C ) = — / dE^ - E, 
m h JEr 



dY{ 



dEp 



v tb v ts 



vr 1 7 1 T S (E C ) 



+ 0{l/m b 2 



(51) 



(52) 



Comparing our result ()52j) for the V u b determination with the one obtained in [HJ [7j one 
expects substantially smaller corrections in the case of Ab baryons, since the corrections 
are only of order \ jm 2 - However, this requires to measure the inclusive semileptonic and 
radiative rare A& decays, which is more difficult experimentally. 
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